Periodic spin textures in a degenerate F — 1 87 Rb spinor Bose gas 
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We report on the spin textures produced by cooling unmagnetized 87 Rb F = 1 spinor gases into 
the regime of quantum degeneracy. At low temperatures, magnetized textures form that break 
translational symmetry and display short-range periodic magnetic order characterized by one- or 
two-dimensional spatial modulations with wavelengths much smaller than the extent of the quasi- 
two-dimensional degenerate gas. Spin textures produced upon cooling spin mixtures with a non-zero 
initial magnetic quadrupole moment also show ferromagnetic order that, at low temperature, coexists 
with the spatially modulated structure. 



Coherent quantum fluids exhibiting spontaneous spa- 
tial order have garnered widespread attention in connec- 
tion to possible supersolid phases of matter pQ and the 
ground states of high-T c superconductors [2] [3] and other 
correlated electronic materials [4 j. Such intrinsically het- 
erogenous quantum fluids may arise due to the interplay 
between multiple order parameters [4j, the influence of 
adjacent ground states with differing tendencies [5] or 
the presence of competing interactions [6j [7] . 

Recent observations hint at similar phenomenology in 
a magnetic quantum gas, the F = 1 spinor Bose gas of 
87 Rb. Early studies of 87 Rb spinor condensates [8j [9] 
suggested that their magnetic properties were governed 
solely by the spin-dependent contact interaction. This in- 
teraction, with mean- field energy — |c2|n(F) 2 , favors spin 
states with maximum magnetization p~Ql[TT]; here, C2 is 
related to s-wave scattering lengths for interatomic colli- 
sions, n is the number density, and F is the dimensionless 
vector spin. Together with the kinetic energy cost for 
spatial variation of the superfluid vector order parame- 
ter, this local interaction favors a simple, homogeneous, 
ferromagnetically ordered spinor condensate. However, 
recent works [T2HT5] point to the significance of mag- 
netic dipolar interactions in determining magnetization 
dynamics in 87 Rb gases with large spatial extent. This 
interaction is long-ranged and spatially anisotropic, and, 
as in magnetic thin films, may favor inhomogeneous and 
spatially ordered spin textures. How the competition be- 
tween the spatially isotropic contact interaction and the 
anisotropic long-range dipole interaction resolves itself in 
such degenerate spinor gases is the subject of several re- 
cent theoretical investigations [16-20 , and remains an 
open experimental question. 

Here, we address this question by examining the mag- 
netic order in gases produced upon cooling unmagnetized 
thermal spin mixtures into the regime of quantum de- 
generacy. At their lowest temperatures, these quantum 
fluids break translational symmetry to arrive at magne- 
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tized spin textures that, while varying between samples, 
consistently display spatial modulations of similar mor- 
phology and length scale independent of initial conditions 
and equilibration time. 

This work differs from previous studies of spin tex- 
tures of spinor Bose condensates [15j [2TJ [22] (also for 
F = 2 gases [23]) in two crucial respects. First, in 
previous works, modulated spin textures arose from dy- 
namical instabilities of long-range-ordered spin-polarized 
Bose-Einstein condensates. Thus, the different character- 
istics of the magnetization structures observed in those 
works, while providing insight on the magnetic interac- 
tions present in the spinor gas, derived from the specific 
unstable initial states chosen. Indeed, the magnetization 
patterns produced through two different dynamic insta- 
bilities, the first being the spin-mixing instability of a 
paramagnetic condensate quenched across a phase tran- 
sition [2U [22] and the second being an instability of a 
helical spin texture [15], showed markedly different spa- 
tial correlations. As such, the present study is aimed 
at revealing the magnetic phases favored intrinsically by 
the spinor gas. Second, this work presents the first ob- 
servation of spin textures in degenerate Bose gases at 
variable temperature, whereas previous works examined 
structure formation only at near-zero temperature. The 
addition of a significant normal gas component should 
add new means for dissipation and relaxation from non- 
equilibrium magnetization states. 

I. EXPERIMENTAL METHOD 

The local one-body density matrix of a spin-1 atomic 
gas describes both rank-1 and rank-2 polarization mo- 
ments, corresponding to the vector spin and magnetic 
quadrupole moments of the atoms, respectively. For 
our experiments, we prepared non-degenerate optically 
trapped gases characterized by homogeneous fractional 
populations (Ci ? Co? C— l) m the three eigenstates of F z . 
These gases were unmagnetized, i.e. characterized by 
zero vector spin, given the dual constraints of £i = C-i 
(zero longitudinal magnetization) and the absence of co- 
herence among the sublevels (zero transverse magneti- 
zation). These constraints still allow for a non-zero 
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rank-2 (magnetic quadrupole) polarization moment, for 
77 = Co — Ci 7^ which breaks spin-rotational symmetry 
by favoring the z axis. 

To prepare such spin mixtures, we began with non- 
degenerate, longitudinally spin-polarized gases trapped 
at the focus of an elliptically focused, linearly polarized, 
825-nm wavelength laser beam. We then produced inco- 
herent spin mixtures by applying resonant tt/2 rf pulses 
while also applying a 50 mG/cm magnetic field gradient 
to the gas. Diffusion of the non-degenerate atoms in the 
inhomogeneous field eliminated transverse coherences, as 
ascertained by probing for Larmor precession in the ther- 
mal gas [24] . Applying a single rf pulse resulted with an 
77 = 1/4 spin mixture, while applying the pulse-diffusion 
sequence repeatedly yielded a fully unpolarized gas with 
77 = 0. 

Following their preparation, spin mixtures were evap- 
oratively cooled by gradually lowering the intensity of 
the optical trapping beam, typically over 200 ms. Dur- 
ing this process, the trap frequencies, with values of 
(uu x , uo y , u z ) = 27r(84, 1000, 10) s _1 at the initial trap 
depth of U/ks — 10 /iK, decreased Vu 
and the gas temperature was found to scale as /c^T = 
0.11 (U/ks — 0.8/iK), with all Zeeman sublevels attain- 
ing the same temperature. The offset in the formula for 
the temperature T accounts for the effects of gravity in 
the y direction. 

After the optical trap depth was reduced to the desired 
level, the spinor gas was allowed to equilibrate, typically 
for another 200 ms, and was then probed by either of 
two methods. The instantaneous vector magnetization 
M = /inF, column-integrated along the y imaging axis, 
was measured in the x-z plane by magnetization-sensitive 
phase contrast imaging, at a measured imaging resolution 
of about 2 fim [24j [25]. Here, ji is the atomic magnetic 
moment and n is the local column density of the gas. Al- 
ternately, to measure the temperature, atom number and 
condensate fraction for each of the Zeeman sublevels, the 
gas was released from the optical trap and subjected to a 
magnetic field gradient that served to separate spatially 
the three Zeeman states. The three components were 
then imaged following an additional time of flight (Fig. 

u 

Throughout the evaporation and equilibration, a B = 
150 mG magnetic field, varying by less than 5 /iG across 
the extent of the gas, was applied along z. This field 
produces a quadratic Zeeman shift of the form qF 2 with 
q/h = (70Hz/G 2 )£ 2 = 1.5 Hz. This shift is smaller 
than the spin-dependent contact and dipolar interactions 
in the degenerate gases studied, with typical energies of 
\c2\n = h x 7 Hz and fiofi 2 n = h x 8 Hz, respectively, 
at a typical density of n = 2 x 10 14 cm -3 . The field also 
induces rapid (110 kHz) Larmor precession of the atomic 
spins, owing to which the magnetic dipole interactions 
assume a precession- averaged form [T4 | \TE \ H7]. 
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FIG. 1. The number of atoms in the thermal fraction of the 
\m z = 0) sublevel (green circles) is compared with the aver- 
age number in the thermal fractions of the \m z = zbl) states 
(blue squares) at different final temperatures of the 77 = 1/4 
clouds. These data are derived from bimodal fits to time-of- 
flight density images taken after the three spin states were 
spatially separated in a magnetic field gradient (such an im- 
age is shown at right). While the relative populations in the 
Zeeman levels at high temperature are consistent with the 
initial value of 77 = 1/4, below T C) jv , the thermal components 
in all Zeeman components become roughly equal. A record of 
the total atom number (condensed and non-condensed) in all 
three sublevels vs. temperature is shown in the inset. 



II. MAGNETIZATION TEXTURES PRODUCED 
FROM 77 = GASES 

Spinor gases produced initially with 77 = correspond 
to the equilibrium state far from quantum degeneracy for 
our system where /c^T ^> {q, |c2|ra}. The bulk features 
of gases produced by cooling this initial mixture display 
several hallmarks of Bose condensation. Gases probed 
after a time of flight show the transition from a Bose- 
enhanced gaussian distribution to a bimodal density dis- 
tribution at a temperature consistent with the ideal-gas 
Bose-Einstein condensation temperature, T c>j /v > given as 
^B^c,Ar = ^(iVo/1.21) 1 / 3 where uj = (cj^c^c^) 1 / 3 and 
No is the atom number in the \m z — 0) state. The 
measured populations within the central peak and the 
gaussian distribution, associated typically with the con- 
densate and thermal fractions, respectively, match closely 
with those expected for scalar Bose gases. The total pop- 
ulation in each Zeeman sublevel remained roughly equal 
as the gas was cooled. 

However, in situ probing of these gases (Fig. [2| reveals 
the formation of complex magnetic structure. Immedi- 
ately below T C) at , the spinor gas becomes spontaneously 
magnetized, breaking spin-rotational symmetry within a 
central region consistent with the spatial dimension of a 
Bose condensate. The magnetization observed for two 
samples produced under similar conditions is shown in 
Fig. [3] These inhomogeneous textures are dominated by 
strong modulation of the magnetization with a charac- 
teristic domain size I of about 5 fim. We observed spin 
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domains of similar size even in gases with reduced atom 
number for which the condensate radii were roughly half 
those shown in Fig. [2] That this length I is smaller than 
the condensate dimensions in the imaged x-z plane and is 
not strongly dependent on the x dimension of the conden- 
sate suggests that the observed translational symmetry 
breaking is an intrinsic tendency of this quantum gas. 

We note also that both I and also the spin healing 
length defined as £ s = h/ (2m\c2\n) 1 ^ 2 = 3 /xm exceed the 
Thomas-Fermi condensate radius along the direction of 
tightest confinement, r y c± 1.5 fim. Thus, spin textures of 
the degenerate gas may be considered two-dimensional. 
In contrast, the scalar healing length, defined convention- 
ally as £ = h/(2mcon) 1 ^ 2 = 0.2 /xm with Co defining the 
spin-independent contact interaction strength, is smaller 
than the Thomas-Fermi radii of the condensate. Thus, 
in terms of scalar excitations, the degenerate gas may be 
considered three-dimensional. 

Importantly, these textures display spatial order, 
characterized by strong modulation along two distinct 
wave vectors. To illustrate this spatial order, we consider 
the spatial correlation function of the vector magnetiza- 
tion, defined as 
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and shown in Fig. [3j The characteristic spatial pattern 
is indicated by the lobes of positive and negative spin 
correlation surrounding the central region of positive cor- 
relation at Sr = 0. This spatial organization is equiva- 
lently indicated by peaks in the Fourier power spectrum 
of the vector magnetization, \M.(k x , k z )\ 2 , at wavevectors 
k ~ tt/1 (Fig.[4| [26]. We highlight three characteristics of 
this spatial order. First, the alternating domain pattern 
is short-ranged, extending regularly only over a typical 
range of 30 /xm that is much shorter than the maximum 
extent (~ 300 /xm at low temperatures) of the magne- 
tized degenerate gas. Second, the modulation pattern 
varies visibly between samples produced at similar tem- 
peratures and equilibration times, with the spin texture 
in some regions characterized by a checkerboard, or two- 
wavevector pattern of magnetization, while in other re- 
gions the modulation along one wavevector appears more 
dominant. In spite of these variations, the orientation of 
the two modulation wavevectors is fairly constant across 
samples. We observe that a rotation of the optical trap 
about y causes a corresponding rotation in the crystal 
wavevectors, demonstrating that the orientation of the 
magneto-crystalline pattern in the x-z plane is pinned 
to the boundaries imposed by the trap. We cannot ac- 
count for a misalignment of these wavevectors from the 
apparent trap axes, though we suspect it arises from trap 
asymmetries due to aberrations in the optical trapping 
laser beam. 

Third, the magnetization modulation of these spin tex- 
tures is characterized by an axis in spin space along which 
the magnetization exhibits the largest variance. To iden- 
tify this local spin axis, we considered the distribution of 
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FIG. 2. Spin textures of 87 Rb F = 1 spinor gases produced by 
cooling thermal spin mixtures with initial values (a) rj = and 
(b) r\ — 1/4. The transverse magnetization M x , y (color wheel 
shown) is shown in the x-z plane. Temperatures are scaled by 
the Bose condensation temperature given the trap parameters 
and the population of the \m z = 0) state. Textures labeled 
(i), (ii) and (iii) are analyzed further in Fig. [4] 



vector spins F = (F x1 F yi F z ) measured at each imaged 
pixel within 30 x 20 /xm 2 regions of the gas (Fig.j5|. These 
regions were centered on the condensate in the x direc- 
tion while the location along z was allowed to vary. The 
observed spin-space distribution was significantly prolate 
[27] . The local spin axis was then defined as lying along 
the largest-eigenvalue eigenvector of the covariance ma- 
trix Dij = ((Fi - (Fi))(Fj - (Fj))) with {i,j} e {x, y, z}. 
This spin axis varies over characteristic distances of ~ 50 
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FIG. 3. The magnetization vector components M XjVjZ 
(grayscale spanning ±2/3 of the maximum condensate magne- 
tization) are shown for the central 25 x 100 /im 2 regions of two 
samples produced by cooling initial mixtures with 77 = to 
80 nK. Regions of sample (a) show prominent magnetization 
modulations along either a single near-vertical wavevector 
(top portion) or along two orthogonal wavevectors (bottom 
portion). In contrast, modulations in (b) are predominantly 
along a near-horizonal wavevector. The magnetization spa- 
tial correlation function G(Sr) evaluated over a 40 x 200 /im 2 
region reveals patterns of alternating spins characteristic of 
these samples. 



fim across the length of the spin texture (Fig. |6|. 



III. MAGNETIZATION TEXTURES 
PRODUCED FROM 77 / GASES 

We also studied the evolution of spin mixtures pre- 
pared with a non-zero initial quadrupole moment (77 = 
1/4). While the small rate of spin-mixing collisions in 
the non-degenerate gas prevented its equilibration to 
rj = within experimentally accessible timescales (sec- 
onds), below T C) at , the thermal fraction of degenerate 
gases did reach roughly equal populations in the three 
Zeeman sublevels (Fig. [I]), presumably due to the high 
density and bosonic enhancement intrinsic to a conden- 
sate. These thermal populations were all consistent with 
the expected non-condensate population for a quantum- 
degenerate, thermally equilibrated, harmonically con- 
fined Bose gas. 

Yet, in spite of this condensate- mediated spin mixing, 
effects of the initial condition were still visible in the mag- 
netization textures of the degenerate gas. In contrast to 
the case of 77 = 0, here we observe two distinct magne- 
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FIG. 4. Spatial Fourier power spectra of the magnetization, 
|M(fcc, k y )\ 2 for gases prepared with 77 = (i) and 77 = 1/4 (ii, 
iii). Periodic spin modulation in (i) and (ii) is indicated by 
broad spectral peaks with modulation wavelength A ~ 10 /mi. 
Ferromagnetic order in (ii) and (iii) is indicated by the spec- 
tral peak at zero wavevector. (ii) and (iii) correspond to gases 
cooled to temperatures of 100 nK and 130 nK respectively, in- 
dicating the abrupt transition to a modulated spin texture. 
The power spectrum is evaluated for the entire magnetization 
texture, corresponding to images shown in Fig. [2] Dashed 
ellipses show demarcation of Fourier space into the central 
region, indicative of ferromagnetic order, and an annular re- 
gion, indicative of short-range spin modulation. Such regions 
are used for an empirical characterization (Fig. |7| of the spin 
textures. 



tization patterns. Just below the condensation tempera- 
ture, the spin textures showed ferromagnetic order, with 
a nearly uniform magnetization that varied only at long 
length-scales of around 100 /im. The stability of these 
transverse ferromagnetic textures shows that the pres- 
ence of a condensate fraction in the \m z = ±1) Zeeman 
sublevels states does not on its own guarantee the appear- 
ance of small spin domains. Only at a distinctly lower 
temperature does the magnetization display periodic spa- 
tial modulations. This modulation pattern coexists with 
ferromagnetic order, constituting a spatial variation of 




FIG. 5. The periodic magnetization textures are charac- 
terized by a spin-space axis along which the magnetiza- 
tion is modulated. Distributions of spins within the central 
30 x 20 /im 2 regions of the magnetization texture, for samples 
at the lowest temperatures. For 77 = 0, both transverse (F x , y ) 
and longitudinal (F z ) spin variations are evident, while for 
77 = 1/4 the modulation is solely transverse. The long axis of 
the prolate spin distribution defines the local spin axis (indi- 
cated by double-headed arrow). For 77 = 1/4, this distribution 
is offset due to the coexistent ferromagnetic order. Spatial 
patterns corresponding to the observed spin-space modula- 
tion are shown schematically at bottom. Dashed ellipses are 
guides to the eye. 
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FIG. 6. The local spin axis, determined with a 30 fim region 
centered variably along the z axis, is not uniform across the 
extent of the condensate but meanders gradually over length 
scales of around 50 /im. The azimuthal and polar angles of 
this axis, given in radians, with typical statistical error bar 
shown, are presented for the lowest-temperature sample in 
Fig.^, with initial rj = 1/4. 



the magnetization about a non-zero average value. 

The distinction between long length-scale ferromag- 
netic order and the short length-scale spin modulation 
is clearly seen in the spatial Fourier power spectra of the 
magnetization (Fig. [4|. This distinction allows us to de- 
fine empirical measures of each type of spatial order. For 
this, we demarcate spatial Fourier space into a central 
region, which we may associate with ferromagnetic or- 
der, surrounded by an annular region which contains the 
broad spectral peaks that characterize the short length- 
scale spin modulation. Summing the spectral power in 
each of these regions provides a measure of the different 
magnetization patterns exhibited by the quantum gas. 

These measures distinguish between spin textures ob- 
tained from differing initial conditions (Fig. [7]). For gases 
produced with 77 = 0, the modulated spin pattern arises 
immediately below T Cj n and intensifies monotonically 
with decreasing temperature, remaining a constant frac- 
tion of the spectral power expected for a fully magnetized 
condensate. The ferromagnetic order remains small at all 
temperatures indicating that the average vector magne- 
tization over length scales larger than ~ 20 /im is roughly 
zero. In contrast, for 77 = 1/4, ferromagnetic order dom- 
inates down to 0.75 T c at = 130 nK. Below this tem- 
perature, the short length-scale spin modulation emerges 
abruptly, as indicated by the different Fourier spectra 
shown in Fig. |4] (ii) and (iii), and coexists with ferromag- 
netism. Here we note that the magnetic dipole energy 
between neighboring spin domains may be approximated 
as fiofi 2 h 2 l/4:7r = x 210 nK. Thus, dipolar interac- 
tions may play a role in the transition toward spatially 
modulated spin textures. 

In addition to the presence of the two distinct mag- 
netization patterns, the textures produced from initial 
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FIG. 7. Measures for ferromagnetic order (red rectangles) 
and for the short length-scale spin modulation (green circles) 
are defined as the sums Pm — |M(fec, k y )\ 2 in regions of 
spatial Fourier space defined in Fig. [4] Insets show these or- 
der parameters as fractions of the integrated spectral power 
expected for a uniformly magnetized condensate, (a) Gases 
cooled from 77 = show a single transition to a magnetic 
texture dominated by short-range spin modulation, (b) Data 
for 77 = 1/4 show an additional transition between textures 
with coexistent ferromagnetic order and spin modulation ob- 
served at low temperature, and solely ferromagnetic textures 
observed at high temperature. 



conditions 77 — and 77 = 1/4 also differ in the overall 
distribution of magnetization orientations. For 77 = 0, 
the magnetization orientation shows no preferred axis or 
plane, while for 77 = 1/4 the magnetization was found to 
lie predominantly in the (transverse) x-y plane. 

This spin-space asymmetry may be explained in part 
by the fact that, for 77 > 0, Bose condensation occurs 
first in the dominant \m z = 0) state. Transverse mag- 
netization may then arise from the dynamical instability 
of this polar-state condensate [21 . This timeline for the 
formation of magnetization was verified by probing gases 
at different times during their evaporation and equilibra- 
tion. Indeed, with 77 = 1/4, unmagnetized condensates in 
the \m z = 0) state were observed before the appearance 
of magnetization. Further, at q = h x 70 Hz > 2|c2|n, 
for which the quadratic shift dominates over both the 
spin-dependent collisional and dipolar interactions and, 
thus, a \m z = 0) condensate is stable, a similar cooling 
sequence resulted neither in a condensate fraction in the 
\m z — ±1) states nor in the appearance of transverse 
magnetization. 
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IV. IMPLICATIONS FOR EQUILIBRIUM 
PHASES 

The limited lifetime (about 1 s) of our gaseous sam- 
ples restricted the duration over which they were allowed 
to equilibrate. This restriction raises the question as to 
what extent the observed magnetic textures display the 
characteristics of thermal equilibrium phases. To help 
address this question, we have verified that both the 
temperature and the populations of the non-condensate 
fractions reach equilibrium within 100 - 150 ms of the es- 
tablishment of the final trap depth. Doubling the evap- 
oration time or lengthening the equilibration period to 
as long as 800 ms produced no discernible changes in 
the magnetic order, other than overall loss in atom num- 
ber. Nevertheless, as discussed above, we observe varia- 
tions in the magnetization patterns in samples produced 
under similar conditions, and persistent differences be- 
tween samples prepared from differing initial spin mix- 
tures. These variations suggest that the magnetization 
of these degenerate spinor gases equilibrates only slowly, 
a property observed also in classical dipolar systems [28] . 
Therefore, the observed textures and the empirical mea- 
sures used to characterize these textures cannot be con- 
sidered to define the equilibrium phase of the degenerate 
spinor gas. 



However, the presence of short-range periodic spin 
modulation is a robust property of the observed spin 
textures, appearing consistently under a variety of ini- 
tial conditions. We have additionally observed similar 
modulation patterns following the dissolution of helical 
spin textures [15] as well as after long evolution times 
following a quench of an unmagnetized condensate to a 
magnetic state (at short evolution times, this instabil- 
ity produces a distinctly different modulation pattern) 
[2lJ[22]. This robustness, and the fact that the character- 
istic length-scale of the spin modulation is constant and 
much smaller than the extent of the gaseous sample (cer- 
tainly in the z direction), supports the claim that such 
structure is a characteristic of the low-temperature equi- 
librium configuration of the bulk quasi-two-dimensional 
spinor Bose gas. 
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